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Abstract 

The NA50 Collaboration has recently observed a strong suppression of J/^f 
production in Pb-Pb collisions at 158 Gev/n. We show that this recent obser- 
vation finds a quantitative explanation in a model which relates the suppres- 
sion mechanism to the local energy density, whose value is higher in Pb-Pb 
collisions than in any other system studied previously. The sensitivity of the 
phenomenon to small changes in the energy density could be suggestive of 
quark-gluon plasma formation. 
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The NA50 Collaboration has recently reported the observation of a strong suppression of 
j/\[r production in Pb-Pb collisions at 158 GeV per nucleon which is not explained by 
conventional models of nuclear absorption. Since such models have been found to account 
reasonably well for all the previous data involving lighter nuclei , the immediate implica- 
tion seems to be that new physics is involved in Pb-Pb collisions, possibly the formation of 
a quark-gluon plasma [|J. 

In this note, we present an interpretation of the data based on the observation that 
the local energy density is higher in Pb-Pb collisions than in any of the systems studied 
previously, in particular the S-U system. In order to explore the possibility that the large 
suppression which is observed in Pb-Pb collisions could be due to the formation of a quark- 
gluon plasma, we adopt a simplified description [[|, in which one considers that all the J/\l/ 's 
produced in a region where the energy density exceeds some critical value are suppressed. We 
emphasize that we do not, and cannot at this stage, make precise statements about detailed 
microscopic mechanisms (for a recent review, see ||). Our purpose is only to test the 
idea that the suppression depends solely on the local energy density, and to see whether 
consequences of this assumption are supported by the data. 

We first review briefly the conventional treatment of nuclear absorption The ratio of 
the J/^/ production cross section in proton- nucleus collisions (ct p a) to that in proton- proton 
collisions (cr pp ) is given by 

M A = i^£d = J_ f d 2 6 (1 - exp (-a a T A (b))) , (1) 

A O'pp AO~ a J 

where Ta(s) = /_ ~ p A (s, z)dz is the nucleon density per unit area in the transverse plane (i.e. 
the plane transverse to the collision axis), and a a is an absorption cross section (although 
we do not write it explicitly, in our calculations a a is multiplied by the correction factor 
(1 — 1 /A), where A is the mass number of the nucleus 0). The quantity Ma m ay be 
interpreted as the probability that a produced J/\l/ survives nuclear absorption. 

In a nucleus-nucleus collision, the survival probability Nab takes, after integration over 
the impact parameter, the factorized form Nab = MaMb- Within the range of A values 
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considered, and for the chosen parametrization of the density, IhNab ~ A 1 ^ + B 1 ^ (see 
Fig.l). We use for p(r) the expression: p^(r)/p = 1/ + exp ( Iz ^ A )) with = 1.1A 1 / 3 
fm, a = 0.53 fm, and po is fixed by the normalization / pA{r)d 3 r = A (e.g. in 208 Pb, 
po = 0.17fm- 3 ). 

As seen on Fig.l, nuclear absorption explains both the proton-nucleus and the nucleus- 
nucleus data up to the S-U system, with a common value of the absorption cross section 
(we have adopted the value a a ~ 6.2 mb used by the NA50 Collaboration [I]]). However, 
the Pb-Pb system deviates significantly from that common trend. One can measure this 
deviation by the ratio 

Nab (measured) 
Nab (estimated) 

where A/ab (estimated) is the value of the survival probability to nuclear absorption alone 
(Nab (estimated) ~ 0.43). The value of as read from Fig.l, is 0.68 ± 0.06. The sur- 
vival probabilities plotted in Fig.l are extracted from absolute cross sections which contain 
systematic errors of the order of 10 to 20%. These errors have been reported in the figure 
although they largely cancel in the relative values of the survival probabilities corresponding 
to a given set of data. In its analysis, the NA50 Collaboration uses the ratio of the J/\l/ pro- 
duction cross section to the Drell-Yan cross section, which reduces most of these systematic 
errors. It also relies, in extracting nuclear absorption, on the transverse energy dependence 
of the S-U data, rather than on integrated data alone. It obtains thus, for the ratio the 
more precise value 0, 72 ± 0.03. 

We turn now to the dependence of the effect on the impact parameter of the collision. 
We write the J/\l/ production cross section at impact parameter b as follows: 

— ^ = T AB (b)N(b). (3) 

where Tab{^) = J d 2 sT^(s) Tg(s — b) is proportional to the probability to produce a cc pair, 
and N(b) is the survival probability at impact parameter b. The quantity N(b) is related 
to Nab introduced above by Nab = (X/AB) J d 2 bTAB(b)N(b). If nuclear absorption is the 
only suppression mechanism, 



m = KUb) = TfAjjr [ d 2 s- 2 (l - e~^«) (l - e -- T ^- b )) . (4) 
1ab{o) J K 

In analogy with eq.(§), we define: 

From the NA50 data, assuming that the last bin in transverse energy corresponds to central 
collisions, i.e. to 6 ~ 0, one extracts the value r^{b = 0) ~ 0.50. 

We now show that the values of these two ratios, defined in eq.(^), and r#(0) defined 
in eq.(||D, can be understood quantitatively if one assumes that the suppression mechanism is 
sensitive only to the local energy density. A central assumption here is that the suppression 
of the J/^f takes place at times short compared with the transverse size of the interaction 
region, i.e. before a substantial transverse expansion of the produced matter has occurred, 
and before the J/^f has traveled a long distance in the transverse direction (we consider the 
production of J/\l/'s near central rapidity). We shall therefore ignore both the transverse 
expansion and the transverse motion of the J/\l/'s. Under these conditions, the fate of a 
J/^f is determined by the properties of the medium in the region where it is created, and is 
controlled by the energy density in the transverse plane. 

To estimate this density, we assume that it is proportional to the density of participants. 
This assumption is motivated by the fact that in nucleus-nucleus collisions, the multiplicity 
and the transverse energy grow approximately linearly with the number of participants ||. 
The participants are the nucleons which collide at least once during the collision of nucleus 
A on nucleus B at impact parameter b. They have a density per unit transverse area given 
by 

n p (s, b) = T A (s) [1 - exp (-cr N T B (s - b))] + T B (s - b) [1 - exp (1 - a N T A (s))} , (6) 

where ~ 32 mb is the nucleon-nucleon inelastic cross section. The total number of 
participants at impact parameter b is N p (b) = J d 2 s n p (s, b). 

A plot of n p for the two systems S-U and Pb-Pb is given in Fig. 2. One sees that, 
up to impact parameters of about 8 fm, there are regions in the Pb-Pb system where the 



density exceeds that in central S-U collisions. The transverse energy dEx/dy achieved in 
central collisions is roughly proportional to N p (0). From this one deduces that the average 
energy density produced in central collisions, proportional to N p (0)/R 2 , is approximately 
the same in the S-U and Pb-Pb systems. However, the maximum density achieved in Pb- 
Pb is about 35% larger than in S-U. One may get an estimate of the maximum value of n p 
by using sharp sphere densities. One gets then n™ ax (6) = 2p ^J {Ra + Rb) 2 — b 2 , which, for 
central collisions, is proportional to A 1 / 3 + B 1 ^ 3 . Note that the the J/^f production, being 
proportional to TaTb, occurs dominantly in the regions of largest density. 

Following [|j], we now model the effect of quark-gluon plasma formation by assuming 
that the J/ty produced at point s is completely destroyed whenever the density at that 
point exceeds a critical value. That is, we calculate 

M{b) = jr^j J d 2 s 1(1- e-- T -( g )) (l - e-- T *( s - b )) 6{n c - n p {s)). (7) 

A plot of M{b) as a function of centrality, defined as N p {b) /N p {b = 0), is shown in Fig. 3, for 
various values of the critical density n c . Given the fact that no suppression is observed in 
S-U collisions other than nuclear absorption, the critical density has to be bigger than the 
highest value attained in S-U collisions, i.e. 3.3 fm~ 2 (see Fig. 2). Choosing this particular 
value for n c , one obtains Nab = 0.28. Therefore, r^, = 0.28/0.43 = 0.66, to be compared 
with the value r^, = 0.72 obtained by NA50. Furthermore, for central collisions, we find 
np(0) = 0.17/0.39 = 0.44, to be compared with the value t>(0) « 0.50 of NA50. Thus the 
two main observations of NA50 can be accounted for quantitatively by this simple picture. 

Within the present model, the value = 0.66 is to be looked at as the lowest possible 
since we have considered the most extreme scenario: total suppression above n c and low- 
est possible value of n c . The fact that the resulting is only slightly smaller than the 
experimental one puts very severe constraints on the suppression mechanism. It seems for 
instance difficult to arrive at a value of as small as 0.72 by a mechanism which would 
set in gradually beyond S-U, and/or suppress only a fraction of the \l/'s, as in a scenario 
advocating the suppression of resonances decaying into J/^, such as x f° r example. In order 
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to attain the value r$ ~ 0.7, all resonances have to be suppressed above n c . 

Since the suppression mechanism is very sensitive to the local energy density, the ob- 
tained value of »> will be also sensitive to a number of factors. For example, it is sensitive to 
the parametrization of the nuclear density. Thus, by taking another common parametriza- 
tion (Ra = 1.19A 1//3 — 1.61A -1 / 3 , a = 0.54 fm) which makes the S nucleus smaller, we get 

= 0.72 instead of 0.66. The deformation of the Uranium nucleus could also slightly alter 
the value of r$. We should also mention that in comparing the two systems S-U and Pb-Pb, 
we have neglected the variation of the colliding energy, from 200 GeV to 160 GeV. Such an 
energy shift results in a slight decrease of the multiplicity density (in p-p collisions, this can 



be estimated from [10] to be about 4%). On the other hand, the energy density is likely to 
increase with the number of participants faster than linearly, as was assumed in our calcu- 
lation. These effects should be taken into account in a more complete calculation. It should 
also be stressed that none of the results on which we rely are totally model independent, 
since the ratios rq, are obtained after extraction of nuclear absorption. 

It would obviously be desirable to get confirmations of the present scenario from indepen- 
dent observables. In particular, it is worth recalling that the existence of a threshold effect is 



not a clear-cut theoretical prediction JTTJ. It is therefore crucial to confirm experimentally 
such an effect and to determine the corresponding density; this could be achieved by explor- 
ing collisions with smaller targets. There are also several effects which could be looked for in 
the Pb-Pb data, and which, if observed, would give confidence in the overall picture. These 
include the effect of fluctuations at large transverse energy, and the expected saturation 
with centrality of the ratio r^i/r^, and of the average transverse momentum squared of the 
J/^'s. We now discuss these three points. 

It can be seen in Fig. 3 that is very sensitive to the value of the critical density. One 
can get a simple estimate by considering sharp sphere nuclear densities, and by neglecting 
nuclear absorption. Then r^(0) = (n c /n™ ax ) 4 for n™ ax > n c , in rough agreement with 
the results displayed in Fig. 3. We can write this same ratio in terms of energy densities: 
r^(0) = (e c /e max ) 4 . Thus, e c being fixed, a 10% increase in e max due to fluctuations, leads to 
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a decrease of about 30% in r^(0). One could therefore observe a further noticeable decrease 
of in collisions involving the largest transverse energies. 

It has been observed in the S-U system that the ratio rv / decreases with increasing 
centrality. Because there are evidences that the J/ty and the suffer the same nuclear 
absorption (JEJ, it is natural to attribute the extra suppression to collisions with comovers 
(it is plausible that the loosely bound ty' is more easily destroyed in hadronic collisions than 
the J/ty)- However, in the present scenario, both the J/ty and the \&' are destroyed before the 
comovers have a chance to do anything. As a result, the ratio rv/r* remains approximately 
constant as a function of centrality, as soon as the critical density is reached, i.e. for b < b c . 
This is easily deduced from Eq.(|7|): when b < b c , r$/(6)/r$(6) ~ J\fcom(b c )J\fw(b c )/Afy(b c ), 
and is approximately independent of b, while it should decrease if no plasma is produced. We 
have made a crude estimate, using the model discussed in |4j] of the quantity N comi which is 
the survival probability of the J/\I/ after its interactions with comovers. The values that we 
obtain depend somewhat on parameters. However, the saturation of the ratio r^'(6)/r$(6) 
with increasing centrality is a fairly robust consequence of the model. 

A natural explanation for the variations of the J/\I/ p^-distributions observed in nuclear 
collisions has been given in terms of initial state scatterings [|13||. In this picture, the increase 



of (pj) at impact parameter b is given by (p?) = (Pt)o + C riAB{b) where C is a constant 
whose value can be determined from proton-nucleus data, and 

n AB {b) = — / d 2 sT A (s)T B (s - b) [T A (s) + T B (s - b)] Af(b) (8) 

J-AB\p) J 

is the average density of nucleons seen by a J/ty. The factor Af(b) has been left out in 
previous analysis. However, it is important whenever the suppression is large. In particular, 
it is responsible here for the fact that n AB {b) remains roughly constant when b < b c , while 
in the absence of a plasma, n AB (b) would continue to increase by some 25%. This result, at 
variance with early expectations that a quark gluon plasma would strongly affect the J/ty 
momentum distribution [13], comes from the fact that in the present scenario all the J/^'s 
need to be suppressed, irrespective of their transverse momentum, when b < b c . 



In conclusion, we have explored a scenario in which J/^f production is totally suppressed 
in regions where the energy density exceeds some critical value. Quantitative agreement 
with the present NA50 data is obtained by choosing the critical density slightly greater than 
the density attained in central S-U collisions. The fact that the maximum densities reached 
in the Pb-Pb and S-U systems may differ by no more than 35% suggests a strong sensitivity 
of the suppression mechanism to small changes in the energy density. It makes it difficult to 
interpret the present Pb-Pb data in terms of collisions with comovers (such an interpretation 
is also made difficult by the theoretical arguments developed in ||15|| ). It is therefore tempting 
to speculate that the large increase in the suppression is due to a dramatic change in the 
properties of the produced matter, pointing to the possible production of the quark-gluon 
plasma. We wish to stress however that the picture presented in this letter is very crude, and 
although it appears to account for the bulk features of the present data, many refinements 
need to be worked out, and detailed confrontations with more data need to be made, before 
unambiguous conclusions can be drawn. 
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FIG. 1. The J/\I/ survival probability after absorption through nuclear matter, as a function 
of A 1 / 3 + B 1 / 3 , where A and B are the mass numbers of the colliding objects. The full line (dotted 
line) is the survival probability for the proton-nucleus (nucleus-nucleus) systems, calculated with 
a cross section a a = 6.2mb. The data at 450 GeV are obtained from refs.fl] and [8], those at 
200 GeV from refs.fl] and [7]. In order to obtain the survival probability, each set of data has been 
rescaled by a constant factor so as to obtain the best fit to the theoretical curve. 
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FIG. 2. The density of participants n p (s), for s along the direction of the impact parameter, for 
various values of the impact parameter: b = 0, 2, 4 • • • fm. The origin is at a distance bj (1 + Rb/Ra) 
from the center of nucleus A. left: S-U collision; right: Pb-Pb collision. The horizontal dashed 
line corresponds to the largest density achieved in the S-U system, n p = 3.3fm~ 2 . 
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FIG. 3. The survival probability of a J/^f in Pb-Pb collisions after absorption in nuclear 
matter and dissolution in a quark-gluon plasma (eq. ([?])). For values n c > 4.4 fm -2 , there is no 
suppression beyond nuclear absorption. The three curves showing an effect of the quark-gluon 
plasma correspond to n c =3.7, 3.5 and 3.3 respectively. The corresponding values of the ratio r$ 
(eq.(^)) are respectively 0.82, 0.74 and 0.66. 
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